Depletion effects in smectic phases of hard rod-hard sphere mixtures 
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It is known that when hard spheres are added to a pure system of hard rods the stability of 
the smectic phase may be greatly enhanced, and that this effect can be rationalised in terms of 
depletion forces. In the present paper we first study the effect of orientational order on depletion 
forces in this particular binary system, comparing our results with those obtained adopting the usual 
approximation of considering the rods parallel and their orientations frozen. We consider mixtures 
with rods of different aspect ratios and spheres of different diameters, and we treat them within 
Onsager theory. Our results indicate that depletion effects, and consequently smectic stability, 
decrease significantly as a result of orientational disorder in the smectic phase when compared with 
corresponding data based on the frozen-orientation approximation. These results are discussed in 
terms of the r parameter, which has been proposed as a convenient measure of depletion strength. 
We present closed expressions for r, and show that it is intimately connected with the depletion 
potential. We then analyse the effect of particle geometry by comparing results pertaining to 
systems of parallel rods of different shapes (spherocylinders, cylinders and parallelepipeds). We 
finally provide results based on the Zwanzig approximation of a Fundamental-Measure density- 
functional theory applied to mixtures of parallelepipeds and cubes of different sizes. In this case, we 
show that the r parameter exhibits a linear asymptotic behaviour in the limit of large values of the 
hard-rod aspect ratio, in conformity with Onsager theory, as well as in the limit of large values of 
the ratio of rod breadth to cube side length, d, in contrast to Onsager approximation, which predicts 
r ~ d 3 . Based on both this result and the Percus-Yevick approximation for the direct correlation 
function for a hard sphere binary mixture in the same limit of infinite asymmetry, we speculate that, 
for spherocylinders and spheres, the r parameter should be of order unity as d tends to infinity. 



I. INTRODUCTION 



In recent years experimental mixtures that closely re- 
semble a hard rod-hard sphere system have been studied; 
typically the rods are represented by tobacco mosaic or 
/d-virus particles, while the spheres are represented by- 
polystyrene latex particles or globular proteins [j], 0, Hf. 
The phase diagrams found include the isotropic phase, as 
well as phases with liquid-crystalline symmetry. Among 
these, nematic (N) and smectic (Sm) phases have been 
found. In addition, bulk dcmixing transitions, as well 
as microsegregated phases of various symmetries, have 
been observed 0. Intensive theoretical effort has been 
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devoted to the understanding of these systems. A useful 
concept in this effort is that of a depletion force. Struc- 
tural and thermodynamic stability of systems consisting 
of hard particles can be understood solely in terms of 
entropic effects which, in mixtures, can be reinterpreted 
as an attractive depletion force. Some studies have been 
done recently in an attempt to quantify attractive deple- 
tion interactions between solute particles, of anisotropic 
shape, mediated by solvent particles that can be isotropic 
or possess themselves liquid-crystalline order 0, M, 111 • 

One of the microsegregated phases that has been ob- 
served is the lamellar phase. This phase, having smectic 
symmetry, consists of alternate pure layers of rods and 
spheres. The lamellar phase can be greatly stabilised 
with respect to the corresponding smectic phase in the 
pure-rod fluid, as predicted in Ref . 7] , and shown in 
Refs. [|, [1 [II m Eld- In the lamellar phase, deple- 
tion forces result from the fact that it is entropically more 
favourable for spheres to occupy the interstitial space, 
which creates a large effectively attractive force between 
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adjacent layers of rods; enhanced smectic-phase stability 
ensues. 

Thus far, theoretical studies of the lamellar phase have 
mostly used, among others, the following approxima- 
tions: parallel rodlike particles with frozen orientational 
order, specific particle geometry, Onsager theory |14| and 
neglect of particle flexibility. Cinacchi et al.[15[ stud- 
ied the hard-spherocylinder (HSPC)-hard-sphere (HS) 
mixture using Onsager theory and free particle orienta- 
tions, but presented phase diagrams including isotropic, 
nematic and the lamellar phase, without specifically ad- 
dressing any stability issue. In this paper we lift stepwise 
the first three of these restrictions and assess their impact 
on the enhancement of smectic stability in the mixture 
with respect to the pure rod system. 

Density-functional theory is a convenient theoretical 
tool to study the structure and phase behaviour of soft- 
condensed matter. Onsager theory is a most success- 
ful version of density-functional theory, intended for 
bodies interacting through hard potentials. Some of 
the shortcomings of Onsager theory are improved by 
Fundamental-Measure theory [l6[ (FMT), which intro- 
duces the fundamental particle measures in the theory 
and should therefore be quite general. Our theoretical 
analysis is based on both these density-functional ap- 
proximations. Specifically, we first use the HSPC model 
within Onsager theory to include the orientational de- 
grees of freedom of the rods by means of an orientational 
distribution function and the associated order parame- 
ter. The results will be discussed in terms of a parame- 
ter, r, which is directly related to the depletion strength 
in smectic phases ■ Then, we discuss the effect of 
particle geometry by considering various shapes for the 
hard particles in the frozen-orientation approximation. 
Finally, we consider a FMT approximation for mixtures 
of hard parallelepipeds (HPAR) and cubes (HC), and 
analyse it in the context of the Zwanzig approximation 
[T7[; the latter consists of considering particle orienta- 
tions to be restricted to three mutually perpendicular 
axes. Recent advances on the application of FMT to 
freely rotating particles have focused on anisotropic par- 
ticles with infinitely narrow breadths (needles). 18] . Since 
one of our aims is to investigate how the depletion mech- 
anism changes as the breadth ratio of the particles is 
varied, these recent developments are not appropriate in 
our context, and the use of the Zwanzig approximation 
seems to be justified for lack of a better approach. FMT 
provides a different theoretical viewpoint with respect to 
which the predictions of Onsager theory can be assessed. 
In fact, the main purpose of this part of the work is to 
assess the impact of a proper inclusion in the theory of 
pair correlations on the depletion mechanism in smectic 
phases. In the framework of FMT, we additionally study 
in detail the asymptotic limit of r with respect to dif- 
ferent parameters relating to the aspect ratio of the rod 
and the relative sizes of HPARs and HCs. With a view 
to completing this study, and based on some expressions 
used to calculate the r parameter, which we evaluate for 



hard spheres, we speculate that, for HSPC-HS mixtures, 
r ~ 1 as the ratio of rod breadth to sphere diameter 
tends to infinity. 

In the following section the general procedure used to 
locate the nematic-smectic spinodal line is illustrated; 
for the sake of clarity, this is illustrated in the context 
of Onsager theory for a binary mixture of hard rods and 
HS's. In Section Hm the t parameter, which is used to 
quantify the stability of the smectic phase when HS's 
are added, is defined, and a closed expression, valid in 
the context of Onsager theory and parallel rods, is pre- 
sented. The results are presented in Section IIV1 which 
has three subsections. The first is devoted to the effect 
of free particle orientations, the second addresses the ef- 
fect of particle shape and the third contains the results 
for mixtures of HPAR's and HCs, analysed by means of 
FMT in the Zwanzig approximation. This last subsection 
is in turn divided into two parts. The first is devoted to 
the results obtained by applying Onsager theory to the 
Zwanzig HPAR-HC system, while in the second we de- 
scribe the results from FMT, together with a detailed 
discussion on the asymptotic limits of the r parameter. 
The conclusions are presented in Section [V] The Ap- 
pendices collect expressions for the Fourier transforms 
of the overlap functions for the various particle geome- 
tries explored in the paper, together with a few details 
on the numerical minimisation of the functional in the 
context of the frozen-orientation approximation. An ex- 
act expression for r in terms of the correlation functions, 
valid for a general mixture of freely rotating particles, 
is also included. Explicit expressions for these functions 
in the Zwanzig-Onsager and Zwanzig-FMT approaches 
for HPAR-HC and HS binary mixtures are provided. 
Finally, the asymptotic behaviour of the r parameter, 
within the framework of FMT, and for the HPAR-HC 
mixture, is described. 



II. ONSAGER THEORY FOR HARD 
ROD-HARD SPHERE MIXTURES AND 
STABILITY ANALYSIS 

The first theory we employ to describe hard rod-hard 
sphere mixtures is Onsager theory, first proposed by 
Onsager 1J| for a pure system of hard rods and exten- 
sively used in theoretical approaches to orientational or- 
dering in hard-rod fluids. It can be regarded as a trun- 
cated virial expansion, up to second order in density, of 
the excess part of the free energy. Thus, it contains the 
orientational-dependent second-virial coefficient exactly. 
It is a density-functional theory, since the free energy de- 
pends functionally on the orientational distribution func- 
tion. The extension of Onsager theory to mixtures is 
straightforward; the reader is referred to Ref.[15] for a 
detailed account on its implementation for the present 
model fluids. 

We consider a two-component mixture of hard rods 
and hard spheres labelled with 1 and 2, respectively. The 
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density functional in the Onsager approximation, assum- The spinodal line is then obtained by solving the equa- 
ing uniaxial smectic symmetry and taking z as the coor- tions 

dinate along the smectic layer normal, is written as r _ -i r 

A(k 0l po,Q )= 1 + poXifn(k ; Q ) 1 + ^0^2/22(^0) 
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where $ = (3J- '/V is the Helmholtz free-energy density 
per unit thermal energy, d s the smectic layer spacing, 
S T (Q) is the orientational-entropy density of the hard 
rods, Q their orientational order parameter, pi(z) and 
pi{z) are local densities for the two components, and 
fij (z) are angular averages of the overlap functions be- 
tween species i and j, already integrated in the xy plane. 
The procedure used to calculate these averages has been 
described in detail elsewhere for the case of mixtures of 
spherocylinders [l5|. In the case where one adopts the 
popular approximation of considering that particles pos- 
sess perfect orientational order (Q = 1), an approxima- 
tion used in all previous theoretical analyses of the hard- 
rod-hard-sphere mixture @, H, H, OH Qjl] , the functions 
fij (z) can usually be written exactly, depending on the 
particle geometry. Note that fu and /12 both depend on 
Q, but certainly not /22- In the nematic phase, pi(z) is 
a constant. 

Since we are interested in searching for the nematic- 
smectic spinodal, we consider the following perturbations 
on the constant nematic densities: 



pi[z) = pXi (1 + Xi cos kz) . 



i= 1,2 



(2) 



where p is the total mean density, x 2 = x, x\ — 1 — x are 
the molar fractions of the two components, and k is the 
smectic wave number. Note that Xi can be positive or 
negative; in the case that A1A2 < 0, smectic layers rich 
in one of the components will alternately appear, giving 
rise to a microsegregated smectic phase. The free energy 
difference between smectic and nematic states, to second 
order in the A^'s, is: 
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where fij(k) are the following cosine Fourier transforms: 
dt/ y (t) cos kt, i,j = 1,2. (4) 

-00 

The Hessian matrix, d 2 ((3AF/N) / dXidXj , necessary to 
study the stability of the nematic phase against smectic 
fluctuations, is 
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for po, kg and Qq, the values of the density, wave- vector 
and order parameter of the unstable nematic at the spin- 
odal, respectively; note that these quantities will depend 
on the value of the composition x. The second equation 
ensures that instability will occur at some particular k 
vector for the first time, as the density is decreased from 
above down to the value at the spinodal. For the pure 
case (x = 0) the instability equations reduce to 



1 + Pohi(ko;Qo) = 0, 
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These are the equations defining a spinodal to an ordered 
phase in a pure system: 1 — poc(fco) = 0, where c(fc) is 
the Fourier transform of the direct correlation function 
(in Onsager theory c(k) = —f(k)). 



III. THE r PARAMETER 

The stability of the nematic phase against smectic-type 
fluctuations can be quantified in different ways. In this 
paper we use the parameter t proposed by Dogic et al. [To| 
in their analysis of the hard-rod-hard-sphere mixture: 



drj 



lim 

V2^o dr]2 



(9) 



In this expression rj is the total packing fraction of the 
mixture, r\i the partial packing fraction of spheres, and 
the derivative is evaluated at the nematic-smectic spin- 
odal line. Other definitions are possible; for example, one 
could use the spinodal line p(rj2), with p the pressure, in- 
stead of 77(772)- Here we adhere to the r parameter as a 
measure of smectic stability. It turns out, as shown be- 
low, that the r parameter can be directly related to the 
depletion potential, so that r contains basic information 
on depletion forces in the system. 

We now obtain a compact expression for r in the frame- 
work of Onsager theory, considering the case of perfect 
orientational order (Q = 1). The generalization to free 
orientational order and a general density functional is 
given in the Appendix C. Our strategy consists of ob- 
taining the nematic-smectic spinodal line 77 = 77(772) per- 
turbatively in the packing fraction r] 2 (i.e. for small val- 
ues of 772), and then extracting the value of r from the 
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first-order term. We begin by first expanding Eqns. (J6j) 
at small x. We assume the expansions 



P = Pa + PqX + 



k = k + k' x + 



(10) 



where now po and are the density and wave-vector, at 
the x = spinodal, and p Q and k' Q are the derivatives of p 
and k with respect to the molar fraction cc, at the x = 
spinodal, respectively. Inserting these equations into the 
first two of Eqns. ([6]) we obtain, at zeroth order in x: 

l + Po/n(fco) = 0, &(fco) = 0. (11) 
To first order in x we get 



Po 



k' = -2^M^l. (12) 
/n(feo)/(i(fco) 



Now, differentiating the equations 7/ = ryi + r/ 2 = p(l — 
x)ui + 1]2 and 772 = pxv2, and taking the limit x — > 0, we 
obtain the relation 
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where vi, V2 are the volumes of the two particle species, 
rods and spheres respectively. This equation, along with 
Eqn. (JIIJ, gives 770 = -Wi//u(fco), and 
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In fact, this expression remains valid for binary mixtures 
composed of any convex bodies, but can only be used in 
the approximation of frozen particle orientations (Q = 1) 
and in the framework of the Onsagcr approximation. The 
r parameter is calculated by solving Eqns. (fTTj) and then 
evaluating Eqn. fT4")) . 

We now proceed to show the direct relation existing 
between the parameter r and the depletion interaction 
between two parallel rods mediated by the solvent parti- 
cles (hard spheres) , at least in the limit where the density 
of the solvent particles p 2 and their diameter, relative to 
the breadth of the rods, are very small. These are the 
conditions under which the Asakura-Oosawa approxima- 
tion [1 91 ] is valid, and the depletion potential becomes 
Vde P (r,P2) = -p2 [fi2 * /12] (r) (the asterisk standing for 
convolution). Now 



and Eqn. 
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(this equation can also be obtained from a mean-field 
perturbative treatment of the fluid, by considering par- 
allel hard rods that interact via a depletion potential in 



mean-field approximation). Eqn. (|16p relates the r co- 
efficient with the depletion potential, and shows that r 
contains basic information on depletion forces. In turn, 
it provides a condition under which r may be negative, 
and thus directly links the depletion potential to the en- 
hancement of smectic-phase stability in fluid mixtures of 
hard rods and spheres. 



IV. RESULTS 

A. Mixtures of HSPC and HS: frozen versus free 
orientations 

All previous theoretical analyses of the lamellar phase 
in the hard-rod-hard-sphere mixtures 0, EH, EH have 
relied on the approximation of considering perfect ori- 
cntational order. In this section we assess the impact 
of this severe approximation on the formation of smec- 
tic phases in the HSPC/HS mixture. We have investi- 
gated a number of HSPC/HS mixtures using the orig- 
inal Onsager theory, changing the HSPC aspect ratio 
n\ = (Li + D\)jD\ and the HSPC breadth relative to 
the HS diameter, d = Dx/Di- Fig. Q] shows the spin- 
odal line 77-772 of the nematic-smectic transition, for the 
cases of frozen and free orientations, and for K\ = 8 and 
various values of the scaled HS diameter. We note that, 
in all cases shown, the smectic phase is stabilised with 
respect to the nematic phase (negative slope at 772 = 0) 
on adding HS to the pure HSPC fluid; this is always the 
case when d > 1 (spheres smaller than the cylinder di- 
ameter). The effect is amplified when both the HSPC 
aspect ratio length (ki) is increased and the sphere di- 
ameter is decreased (d increases). However, the effect is 
less pronounced when the HSPC are free to orient their 
main axes. The r parameter, shown in Fig. reflects 
this behaviour, which points to a less strong depletion 
effect due to orientational fluctuations. 



B. Other hard rod— hard sphere mixtures: effect of 
particle shape 

To investigate the effect of the shape of the rods on the 
depletion effect and, in turn, on the stability of the smec- 
tic phase, Vesely[13| has recently considered parallel- 
rod models consisting of linear overlapping hard-sphere 
chains, hard ellipsoids and hard spheroellipsoids. Here 
we have studied, still in the approximation of frozen ori- 
entations, two additional types of hard particles: cylin- 
ders (HCYL) and parallelepipeds with a square trans- 
verse shape (HPAR). To minimize the effect of the dif- 
ferent particle geometries, the particles were chosen to 
have the same volume (see later). In all cases the minor- 
ity component continues to consist of hard spheres. The 
functions fij(z), and also their Fourier transforms, can 
be calculated analytically in these cases. The relevant 
expressions can be found in the Appendix A. 
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FIG. 1: Nematic-smectic spinodal line 77(772) , according to 
Onsager theory, for the case Ki = 8, and d _1 = 1,1/2 and 1/3 
(indicated as labels), (a) Frozen-orientation approximation; 
(b) free orientations. Note that horizontal and vertical scales 
in both graphs are the same. 



Results for the r parameter in the case of HCYL/HS 
and HPAR/HS mixtures, along with the previously con- 
sidered HSPC/HS mixtures, are shown in Figs. HKa)-(c). 
The t parameter is plotted against the aspect ratio of the 
rods Ki in Fig. [3{a) for different values of the scaled in- 
verse HS diameter d = D\ / D2 (D\ being the side length 
of the parallelepiped in the case of HPAR); Fig. [3Jb) is 
a zoom of (a) in the region of small K\. Finally, the r 
parameter against d for different values «i is plotted in 
Fig. G^c). We can see from Fig. E^a) that r tends to 
decrease linearly, becoming negative for sufficiently large 
aspect ratios. This effect is more and more pronounced 
as the diameter of the HS is decreased. An interesting 
feature is that, when the diameter of HS is greater than 
that of the rod, depletion is larger, and therefore smec- 
tic stability is enhanced more substantially, as one goes 




FIG. 2: Values of r as a function of inverse HS diameter ratio 
d, for various values of HSPC length-to-breadth ratio m. 
Lines: Frozen-orientation approximation (Q = 0); from top 
to bottom: k\ = 6, 8, 11 and 21. Symbols: free orientations. 
Squares: fti = 5; open circles: ki = 8; filled circles: n\ — 11, 
and triangles: Ki = 21. 



from HPAR to HCYL or HSPC (note that these two are 
very similar). The opposite behaviour (i.e. depletion in 
HPAR is enhanced with respect to HCYL and HSPC) 
results when the HS diameter is equal or less than that 
of the HSPC [Fig. G^a)]. For short rods r exhibits oscil- 
latory behaviour [see Fig. E^b)]. Finally, fixing ki and 
increasing d, we find that the r parameter decreases as 
a cubic power-law in d and that its value for HPAR is 
much less than that for HSPC or HCYL. 

The linear and cubic dependences of t with respect to 
K\ and d, respectively, can be elucidated from Eqn. |T5|) 
using the fact that the volume ratio can be expressed as 
V1/V2 — Kid 3 for the HCYL case (the other cases being 
similar), while the squared ratio of Fourier transforms of 
the overlap functions tends asymptotically, for large k\ 
and d, to a constant. Simple calculations give, in this 
limit. 
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= tt: cos -7- ) K\d' 



(17) 



with fcg = koLi the reduced wave number of the N-Sm 
spinodal instability in the one-component system. 

To study the effect of a varying particle shape on the 
phase behaviour of the rod-sphere mixtures, we have per- 
formed a full minimization of the Helmholtz free-energy 
functional with respect to smectic-like density profiles of 
the rods [/Oi(.z)] and spheres [p2(z)]- The minimization 
was carried out with respect to the Fourier amplitudes 
of the truncated Fourier expansion of the density pro- 
files and smectic wave number (details are provided in 
Appendix B). The choice of a Fourier expansion is justi- 
fied by the fact that the HS equilibrium density profiles 
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FIG. 3: (a) The r parameter vs. ki for different values of 
the parameter d (indicated in the graph), (b) A zoom of 
(a) for small /ti's. (c) r vs. d for different values of Ki as 
indicated in the figure, (c) r parameter as a function of d for 
different values of Ki, as indicated. Results for HSPC, HCYL, 
and HPAR are plotted with solid, dashed, and dotted lines, 
respectively. 
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FIG. 4: Equilibrium density profiles of the equimolar HSPC- 
HS mixture (different components are labelled in the figure) 
resulting from minimization with respect to the Fourier am- 
plitudes (solid lines) and with respect to a single-amplitude 
exponential parametrization (dashed lines). The HSPC has 
an aspect ratio ki = 10 and its breadth coincides with the HS 
diameter, taken as unity. The pressure in reduced units was 
chosen to be 0.4. 



Since we would like to compare the results from the 
density-functional minimization of mixtures composed of 
rods of different geometries (HSPC, HCYL and HPAR), 
a sensible criterion is required to choose the sizes of the 
particles so as to make sure that these differences arise 
from the shape and not from different sizes and volumes. 
For this purpose we have used three different criteria: 
(i) all particles have the same diameter and the same 
length, (ii) all particles have the same length and the 
same volume, and (iii) all particles have the same aspect 
ratio and the same volume. Comparing the equilibrium 
density profiles resulting from use of these three criteria, 
we have concluded that differences between particles are 
minimised using the third criterion (these results are not 
shown here). 

We have calculated, using the free, i.e. Fourier-based, 
minimization, the phase diagrams of binary mixtures of 
rods (HSPC, HCYL and HPAR) and HS. The HSPC and 
HS diameters were taken to be equal and adopted as unit 
of length. The aspect ratios of all the rods were set to 
10, all particle volumes being the same. The following 
relations result: 



D H c — 1 (3Khspc) 



-1/3 



Dr, 



obtained with this minimization differ very significantly 
from those obtained via the usual single-amplitude ex- 
ponential parametrization, as shown in Fig. |4j We can 
see that the parametrization produces too sharp density 
peaks in the HS profile, which enhances the depletion 
effect. 



-Dhp = ^-jj D HC , = k hspc D^, (18) 

with /i=HSPC, HC or HP. The resulting phase diagrams 
are shown in Fig. [5] The first interesting feature of 
all phase diagrams is the presence of a tricritical point, 
where the nature of the N-Sm transition changes from 
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FIG. 5: Phase diagrams in the reduced pressure-composition 
plane of three different mixtures composed of rods of different 
geometries (HSPC, HCYL, and HPAR) with HS. Different 
lines corresponding to each mixture are labelled in the figure. 
Also indicated are the stability regions of the nematic and the 
smectic phases. Solid lines represent the binodals of the N- 
Sm coexistence, while the continuous N-Sm transitions have 
been indicated by a dashed line. 



elucidate the effect of restricted orientations separately 
from that inherent to the proper inclusion of higher-order 
correlations, a study of the HPAR-HC mixture using a 
Zwanzig-Onsager theory, together with that on mixtures 
of parallel HPAR and HC, is necessary. 

1. Zwanzig model in the Onsager formulation 

In Appendix D we have written expressions for the 
Fourier transforms of the correlations functions Cy (k) as 
obtained from the Zwanzig-Onsager approach. These 
expressions are necessary to calculate the r parameter 
using Eqn. (jC.lOjl . The results obtained are plotted in 
Fig. O In the same figure the parallel case is also plotted 
for the sake of comparison. The usual cubic power-law 
behaviour of r with respect to d is obtained; this is typical 
of any Onsager approach, as discussed in Sec. IIVB1 The 
Zwanzig approach gives a lower value for r as compared 
to the parallel case. This result is similar to that already 
obtained using the freely rotating Onsager approach for 
HSPC. For high values of K\, both approaches collapse 
into a single line, which is due to the high value of the 
order parameter Q of a fluid composed of long particles. 



second to first order [20]. For pressures higher than that 
of the tricritical point, the demixing instability region 
becomes wider. The width of this demixing window in- 
creases for mixtures going from HSPC to HCYL and 
then HPAR. This scenario clearly demonstrates that in- 
creasing differences in shape between spheres and rods 
enhance the depletion in smectics, resulting in a wider 
demixing gap. 



C. Zwanzig model: Onsager theory versus FMT 

In order to check whether the effects presented above 
are robust with respect to the theory and approxima- 
tions used, we have analysed a mixture of HPAR and HC 
using Fundamental-Measure theory in the Zwanzig ap- 
proximation. Zwanzig (2lj introduced a model that con- 
siderably reduces the complexities associated with min- 
imising a density-functional that depends on particle ori- 
entations. The approximation consists of restricting the 
possible orientations of a rod to lie along the three or- 
thogonal axes x, y and z, treating a pure system as a 
three-component mixture (one component for each orien- 
tation). The approximation is very useful and has been 
used quite substantially. In our FMT context, the re- 
sulting theory treats spatial correlations very accurately 
so that it should better describe the onset of smectic or- 
der in the nematic fluid. We have used the HPAR-HC 
system and the Zwanzig approximation because in this 
combination a FMT can be rigorously formulated [l7j |. 
The Zwanzig approximation should not be terribly lim- 
iting, since the nematic fluid is well oriented. In order to 
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FIG. 6: r parameter as a function of d for three different 
values of n\, as indicated. Solid and dashed lines correspond 
to the Zwanzig and parallel approaches, respectively. 



2. Zwanzig model in the FMT formulation 

Details on this approximation can be found in Ref. 
[l7| . Here we only give a brief sketch of the theory. We 
continue to use the notation introduced in the previous 
section for the dimensions of the particles, i.e. the length 
of the HPAR is L\, the side length is D\ (these particles 
are assumed to have a square section) , and D2 is the side 
length of the cubes. Since the unit vector SI only has 
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three possible orientations, the one-particle distribution 
functions p s (r, ft) can be expressed as 



»(*"> O) = ^2 p SIJ ,(r)5(£l - e M ) 



(19) 



where e M , // 



1,2,3, are unit vectors along the three 
perpendicular directions xyz, respectively, and p S /x(r) are 
the local density of species s parallel to the /i-axis. The 
excess part of the free-energy density, in reduced thermal 
units, is obtained in |17| . and has the form 



$ ox (r;{p s }) = -n ln(l - n 3 ) 



ni • n 2 



n 2x n 2 yn2z 



l-n 3 (l~n 3 y 



(20) 



with the functions {n a } (a — {0, la;, ly, lz, 2x, 2y, 2z, 3}) 
being weighted densities obtained as 



n a (r) 



2 3 

EE 

S — l 



dv'p s ,{v>)J${v 



r% (21) 



where 



,(<*) 



are characteristic functions whose spatial in- 



tegrals give the fundamental measures of the particles 
(edge length, surface and volume). The ideal part of 
the free energy density in reduced thermal units for this 
model is 



(22) 



so that the total free energy per unit volume and unit 
thermal energy can be calculated as 



1 r B 

* = T / dz I $ id(z) + $exc(z) 

d s Jo 



(23) 



Density and order-parameter profiles can be defined; in 
particular, the density profile is 



Ps{z) 



Ps^{z). 



(24) 



The nematic-smectic spinodal line can be obtained in 
this theory using the same kind of arguments used for 
the Onsager theory (Sec. HI|) . Note, however, that the 
expression (I14[) for the r parameter is not valid in this 
context, and we need to use the more general formula 
(IC.10[) . valid for any mixture of freely rotating particles. 
In the Appendix C we describe in detail how this formula 
is obtained. An alternative way to calculate this param- 
eter is to implement a numerical differentiation scheme 
on the spinodal line 77(772). 

The nematic-smectic spinodal lines for the case K\ = 8 
and various values of the scaled cube side length d^ 1 = 
D2/D1 are shown in Fig. [Jj whereas the corresponding r 
parameters for different values of K\ are shown in Fig. [5] 
As was the case in the mixtures analysed using Onsager 
theory, r decreases as the size of the cubes is diminished, 
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FIG. 7: Nematic-smectic spinodal line 77—772 , as obtained from 
Fundamental-Measure theory for the HPAR-HC mixture, for 
Ki — 8 and various values of d~ , as indicated in the graph. 



indicating a stronger depletion effect and a corresponding 
enhancement of the smectic-phase stability. An interest- 
ing feature of these results is that they appear to exhibit 
linear behaviour for large values of d (i.e. as the size of 
the cubes becomes smaller). Fig. [8] shows this by means 
of least-square fits to linear functions. The coefficients of 
the fit are functions of the aspect ratio k 1; and, in turn, 
appear to be linear in K\, as shown in Fig. [9] 
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FIG. 8: t parameter as a function of d for various values of 
the parallelepiped length-to-breadth «i = L1/D1 (indicated 
in the graph) . The data are indicated by symbols; linear least- 
square fits r = q(ki) + /3(fci)d are included as lines for each 
value of K\. 

To explain the linear behaviour obtained for the r pa- 
rameter as a function of d, we have analytically calcu- 
lated the asymptotic behaviour of r in the regime d ^> 1 
for a mixture of parallel parallelepipeds and cubes in the 
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FIG. 9: Symbols: coefficients a (circles) and f3 (squares) as a 
function of ki. Solid lines: linear fits a(/ci) = ai + ct2^i and 
/3(ki) = /3i + /3 2 ki, with #2 = -0.038. 



FMT formulation. In Fig. [Til] we compare the r param- 
eter, for various values of K\, obtained from the Zwanzig 
and parallel FMT approaches. An interesting feature is 
that depletion is enhanced in the Zwanzig approach, a 
result opposite to that from the Onsager approach (Fig. 

El). 
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the two terms in the numerator are of order d~ 2 . Since 
the volume ratio is t>i/i>2 = Kid 3 , we obtain the asymp- 
totic behaviour r ~ Kid = L\jDi. This result can also 
be obtained using the Zwanzig model (Figs. [5] and [5]), 
which shows that the third term of the numerator in Eqn. 
(|C.10|) is also of order ~ d~ 2 . The coefficient /3(ki) of the 
linear fits of r with respect to Ki (Fig. [9]) depends lin- 
early on Ki and has a slope f3^ = —0.038 (see caption of 
Fig. [9]), which is in the range of values predicted by the 
asymptotic limit of W(ki) = r/(Kid) as a function of K\, 
shown in Fig. 

The function (dA/dp 2 ) (0) [Eqn. (jFTT]) ] is of order d~ 2 
only within the FMT formulation, the Onsager approach 
giving a term ~ dr. This, in turn, shows the importance 
of properly taking account of pair correlations between 
particles in order to adequately describe depiction inter- 
actions. The asymptotic limit coincides with the adhe- 
sive limit because, when the mixture is highly asymmet- 
ric, the attractive depletion potential between two big 
particles becomes narrower and deeper, tending in the 
limit to a Dirac delta function at contact. This result 
is confirmed by our calculations in Appendix F, which 
show that the inverse Fourier transform of one of the 
terms of the function (dA/dp2)^ is a Dirac delta func- 
tion at contact. We have repeated the same analysis 
in Appendix F for the case of a HS binary mixture in 
the FMT formulation (PY approach). We have obtained 
that (dA/dp 2 )^ ~ 0(d- 3 ). Also, a Dirac delta func- 
tion at contact is present in the inverse Fourier trans- 
form of (dA/dp2)^ ■ The difference in the limit of in- 
finite size asymmetry between the HPAR-HC and HS 
cases is due to the difference in particle shapes: when 
two big particles are closer enough the parallel sides of 
two parallelepipeds exclude many more small particles 
from their interstitial space than the curved surface of 
two big hard spheres. The discussion above forces us to 
speculate about the possible asymptotic behaviour of the 
r parameter for the HSPC-HS mixture in the limit d 3> 1: 
since vi/v 2 ~ d 3 , r ~ 0(d°) [Eqns. (|CT0ll and (jR2l) ]: 
this is a consequence of the caps of the spherocylindcrs 
being spherical. This argument does not apply to the 
behaviour of r in the limit of large K\ for fixed d, which 
continues to be linear. 



CONCLUSIONS AND PERSPECTIVES 



FIG. 10: The r parameter as a function of d calculated from 
the Zwanzig (solid line) and parallel (dashed line) FMT ap- 
proaches. Values of Ki are indicated as labels. 

Details on the calculation of the asymptotic behaviour 
of t with respect to d have been relegated to Appendix 
F. From the general expression of r (Eqn. IC.lOp . partic- 
ularized for the parallel case (Q = 1, which makes the 
third terms in the numerator and denominator to van- 
ish) , the linear behaviour can be explained if the sum of 



The first result from this study on depletion effects in 
smectic phases is that, within Onsager theory, the in- 
clusion of orientational degrees of freedom significantly 
reduces the enhanced stability of the smectic phase over 
the nematic phase as compared to that obtained from 
the frozen-orientation approximation. A second conclu- 
sion that can be drawn is that when the surface curvature 
of the rods is very different from that of the other species 
(e.g. the case of HPAR and HS), depletion is enhanced 
and the smectic phase stabilizes at lower packing frac- 
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tions as HS become more abundant. This effect is also 
reflected in the fact that the N-Sm demixing gap in the 
phase diagram widens as one goes from HSPC to HP 
particles through the HC geometry. Finally, the inclu- 
sion of more realistic pair correlations between particles, 
in the manner of FMT theory, changes the asymptotic 
behaviour of the r parameter (a convenient measure of 
depletion in smectics) for large values of the ratio d, i.e. 
the asymmetry between the breadth of the rod and the 
sphere diameter; in particular, for HPAR, r increases lin- 
early with d, while for HSPC we speculate that r goes 
to a finite limit. The limiting case of large aspect ratios 
and finite d is similarly captured by both Onsager and 
Fundamental-Measure theories. 

In this work we have studied the depletion mecha- 
nism along the N-Sm spinodal line for different mixtures 
by explicit calculation of the r parameter evaluated on 
this line in the pure fluid. An interesting task could 
be to extend this study to equilibrium smectic phases, 
consisting of well-developed density peaks, by defining 
some quantity evaluated at the smectic density profile. 
In this case depletion-based mechanisms may add inter- 
esting phenomenology, such as strong microsegregation 
between different species, which could be directly quan- 
tified using some suitably defined parameter. Work along 
this direction is currently in progress. 



forms are 

/11(A) = J{sin[A(£)i+Li)] 

- Difccos[fc(Di +Li)] -sin(fcii)}, 

/12(A) = Jf{sin[(£>i 2 +Li)fc/2] 

- ^D 12 cos [(D 12 + Li)fc/2] - sin (fcLi/2) j , 

4-7T 

/22(A) = ^3 [sin (kD 2 ) - kD 2 cos (kD 2 )\ , (A.l) 

where D\ 2 = D\ + D 2 . 

For mixtures of HCYL of length L\ and breadth Di, 
and HS of diameter D 2 , the appropriate Fourier trans- 
forms: 
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VI. APPENDIX 



Appendix A. Overlap functions: case of frozen 
orientations 



'kLA (kD 2 



sm 



/22(A) - — [sin (kD 2 ) - kD 2 cos (kD 2 )} 



kU 



(A.2) 



When particle orientations are parallel, the functions 
fij (z) and their Fourier transforms can be written exactly 
for a large class of particle shapes and their mixtures. We 
consider in turn each of the particle geometries analysed 
in the paper. 

In the case of mixtures of HSPC of length L\ and 
breadth D\, and HS of diameter D 2 , the Fourier trans- 



where J\{x) and H\(x) are the Bessel and Struve func- 
tions of first order, respectively. 

For mixtures of HPAR with length and square side 
length Li, D\, respectively, and HS of diameter D 2 , the 
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Fourier transforms are 



fu(k) 



- sin (kLi) , 



/12(A) = 8 



£>xV sin [k(L 1 +D 2 )/2] 
k 



To calculate phase equilibria in the mixture it is more 
convenient to work in conditions of constant pressure. 
Then, we minimize the Gibbs free energy per particle 
and unit thermal energy, 



P 



(B.5) 



7T 
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Ih 

k 2 

D1D2 
2k 



k 3 



cos 



fc(il+X>2) 



kLi 
sin — 



k(Lx + D 2 ) 



kLA fkD 2 
Jl { — 



'kLA (kD 2 
[ —)"'{ — 



47T 



fcil 



/ 22 (fc) = T3 [sin (kD 2 ) - kD 2 cos (kD 2 )) . 



(A.3) 



Appendix B. Minimisation of Onsager functional in 
the frozen-orientation approximation 

We assume a Fourier expansion for the density profiles 
Pi(z) = pxii>i{z) = pxj I 1 + s£> cosnkz , (B.l) 



with fc the smectic wave number, and s„ the n-th Fourier 
amplitude for species i. Introducing this expression into 
the free-energy funtional, Eqn. ([T]), we obtain 



$ = pi \ogp - 1 + ^x.ln 



F — / (z) lni/'i (z) 



+ pW(o:,fc,{4 l) }) r , 



(B.2) 



Here we have expressed the excess part in terms of the 
function 



W(x.k,{s^}) = i£ 



1 -J 



1 Ctt 



(B.3) 



where k n = nk and the Fourier transforms f%j(k) were 
defined in Section |TT] (since we are assuming perfect ori- 
entational order, we set Q = 1). From this, the pressure 
is 



pP = P + p 2 W(x,k,{sH ) }) 



(B.4) 



with respect to the smectic wave number and the co- 
efficients {sn }• The corresponding derivatives can be 
written down explicitely but they have to be solved nu- 
merically using an iterative method. In practice, the 
Fourier expansions have been truncated; they include ca. 
40 terms so as to satisfy a stringent convergence criterion 
in the iterative procedure (the number of Fourier ampli- 
tudes was chosen to guarantee an absolute error less than 
10~ 7 in the density profiles). 



Appendix C. r parameter: general case 

The strategy is to obtain the nematic smectic spinodal 
line perturbatively in the composition x, since we only 
need to know the spinodal in the neighbourhood of x = 
to obtain the r parameter. The equations to solve are 
three: two defining the spinodal, the third giving the 
equilibrium state of the nematic fluid on the spinodal 
[Eqns. (JBJ)]. These equations can be written as 

A(fc,pi,/52,<2) = [1 - pic n (k,pi,p2,Q)] 

x [1 - P2c 22 (k,pi,p 2 ,Q)} - pip 2 c 2 12 (k,p 1 ,p 2 ,Q) = 



dA(k,p 1 ,p 2 ,Q) 
dk 

d$( Pl ,p 2 ,Q) 



dQ 



0, 



(C.l) 



(C.2) 



(C.3) 



where 5y (fc, pi, p 2 ,Q) are the Fourier transforms of the 
direct correlation functions between species i and j. The 
relationship between infinitesimal changes in the vari- 
ables {k, pi, p 2l Q} along the spinodal can be calculated 
from (jC.lj) as 



dA\ , fdA\ , fdA\ _ (3A\ ,, 



(C.4) 



Similarly, the relation between the changes in the order 
parameter Q and densities {pi, p 2 } can be obtained from 
(JC3| as 



[^J dpi + [mp-J dP2 + {w ]dQ = - 



(C.5) 
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Defining the coefficient v 
Eqns. (JU3) and (TC31) : 



dpi/dp 2 , we obtain from Eqn. (|C.10|) . the r parameter can be re-expressed as 



dA 

dp 2 



8A\ (d 2 1>/dQdp 2 ) 
dQJ (<9 2 $/<9Q 2 ) 



dA 
dpi 



8A\ (d 2 <$>/dQdpi) 1 
dQ J (<9 2 $/<9Q 2 ) 



(C.6) 



where the condition given in (|C.2[) of the absolute mini- 
mum of A(fc, pi, P2, Q) with respect to the wave number 
k was used. Note that this expression is valid for any 
composition x of the mixture. Now, evaluating all the 
derivatives at p 2 = we obtain, from Eqn. (|C.1|) 



dpi) 

9A\ (0) = 
dp 2 J 

dA\ {0) = 
dQJ 



-Pa 



-pa 



( dc 



li 



V dpi 
( 

\ dp2 



(0) 



(0) 



'5cnV 0) 



(C.7) 
(C.8) 
(C.9) 



where the superscripts mean that all derivatives are eval- 
uated at x — 0, k = ko, and Q = Qo, values correspond- 
ing to the one-component fluid spinodal [note that, from 
Eqn. (|C.1|) . we get 1 — poCu = at p 2 = 0, which has 
been used to obtain the derivatives above]. 

Now, using the fact that, for the one-component fluid, 
1 — po&ii — 0, and that the r parameter can be expressed 
through the coefficient v evaluated at i = as t = 1 + 
{vi/v 2 )v (°\ we finally obtain 



Vl 

v 2 



X~ \ (o) 
dc n x 

dp 2 



(0) 



gcn\ w (d 2 §/d P2 dQ) 
dQ J (<9 2 $/<9Q 2 ) (0) 



dpi 



dhiY ' (d 2 $/dpidQ) 
dQ 



(<>) 



(«9 2 $/«9Q 2 ) (0) 
(CIO) 



which is a general expression for the r parameter. In 
the Onsager approximation, we have cy(fc, pi, p 2 ,Q) = 
fij{k,Q) and the derivatives of with respect to the 
densities pi vanish, while for parallel rods (Q — 1) the 
derivative with respect to Q also vanishes, and we obtain 
expression (TH|. 

For a fluid mixture of particles without orientational 
degrees of freedom, such as a mixture of parallel paral- 
lelepipeds and cubes or a mixture of hard spheres, expres- 
sion (|C.10|) for r can be reinterpreted as follows. From 



T = 1 



t'l 



'd_s^_ 

dp 2 



(o) 



V2 ( es-f 

\ dpi 



(o)- 



(C.ll) 



where S~^(k, pi, p 2 ) is the inverse structure factor of an 
effective one-component fluid of particles, labelled as 1, 
with interactions between them being mediated by parti- 
cles labelled as 2. This structure factor can be calculated 
by evaluating the second functional derivative with re- 
spect to Pi(r) of the semi-grand canonical free-energy 
functional 



T[pi] =F[pi,p 2 ] ~ p 2 drp 2 (r), 



(C.12) 



at the bulk densities [22|. Here F[pi, p 2 ] is the Helmholtz 
free-energy functional of the binary mixture and the den- 
sity profile p 2 (r) in Eqn. (|CJ. 12|1 is calculated from the 
condition of fixed chemical potential of species 2: 



SF[pi,p 2 ] 



P2- 



5p 2 (r) 

The result for this effective structure factor is 



(C.13) 



S e f f {k,pi,p 2 ) 



1 - picu(k,pi,p 2 ) 
Pip 2 c\ 2 {k,pi,p 2 ) 
1 - P2C 22 (k,pi,p 2 )' 



(C.14) 



After the inclusion of its first derivative with respect to 
pi and p 2 , evaluated at p 2 = in Eqn. (|C.11|) . we exactly 
obtain Eqn. (|C.10[) without the third terms in both nu- 
merator and denominator (as they vanish for fluids with- 
out orientational degrees of freedom). 



Appendix D. Correlation functions for HPAR-HC 
mixture in the Zwanzig— Onsager approach 

In the Zwanzig approximation the rods point along one 
of the Cartesian axes x,y or z. For the binary HP-HC 
mixture in the Zwanzig approach the correlation func- 
tions evaluated at k = (0, 0, k) can be calculated as 

- Cn(k) = 2x]_ fl x ,lx(k) + flx,ly(k) +x\ifu,u(k) 

+ ixxx\\fix,iz(k), (D.l) 



-5i 2 (fc) = 2xxfix,a(k)+x\\fi Zt 2(k), 



-c 22 (k) = f 22 (k), 



(D.2) 



(D.3) 



where the subindexes Ip (p = x,y,z) label particle 
1 (rods), which point along the direction p, while the 
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subindex 2 labels cubes. x± and xu are respectively the 
fraction of rods perpendicular and parallel to the nematic 
director (which is taken to be parallel to the z axis). 
These variables are functions of the nematic order pa- 
rameter Q [i.e. x± = (1 — Q)/3, and xu = (1 + 2Q)/3], 
the equilibrium values of which should be calculated from 
the extremum condition of the free-energy density with 
respect to Q. This condition reads 



dQ 



:Pi 1 1n 



1 + 2Q 
1-0 



-2 Pl (Li-Di) 'DiQ} = 

(D.4) 



The expressions for the Fourier-transformed overlap func- 
tions are 



hx,iy{k) = 2{L 1 +D 1 



2 sin(fcA) 



sin(fcA) 



/ix,i*(A)=4A(A + A 



sin [k(L! + A)/2] 



/i*,i*0) = 8Ui 



2 sin(fcLi) 



A: 



/i X)2 (A)=2(A+A)(A+ A 



^in [fc(A + D 2 )/2] 



/i*, 2 (A) = 2(A+A) 



2 sin [fe(Li + AO/2] 



(D.5) 



Note that, due to the discrete axial symmetry of the par- 
ticles, some overlap functions can be expressed in terms of 
others; for instance, fi x> \ x = fi y ,i y , and fu,ix = fiz.iy 



Appendix E. Correlation functions for HPAR-HC 
and HS mixtures in the FMT approach 

To obtain the Fourier transforms of the correlation 
functions for HP-HC mixtures and for HS mixtures 
we have used the FMT functional, which has the 
same structure in both systems, namely: 



cy(k) 



- v (o) 



/ y -(k) + x ( 1 )^(k) + x ( 2 )^(k) 



+ X (3) ^(k), 



(E.l) 



where , Aj , Sij and Vij are Fourier transforms of the 
overlap function, and of the mean radius, surface and 
volume of the overlap region between particles i and j, 
respectively. For a HPAR-HC mixture the mean radius 
and the surface area of the overlap region are vectors 
oriented along the directions parallel to the edge lengths 
(mean radius) and perpendicular to the sides (surface 
area) of the parallelepipeds. For the HS mixture they are 
scalars, as are the quantities x and x^ 2 \ which for the 
HPAR-HC mixture are vectors. Finally, the correlation 



functions depend on the wave vector k or on its mod- 
ule k in the HPAR-HC and HS mixtures, respectively. 
The expressions for \ corresponding to the HPAR-HC 
mixture are 
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with p = pi and 77 = piVi the total 
density and packing fraction, respectively, £ 2 — 

i&y&zi&z&xi&x&y), and 

gls = £iy = ^PiA) £iz = P1L1 + p 2 D 2 , (E.5) 

i 

6, = = PiAA + p 2 A!, ^ 2 = X)«A 2 , (E.6) 

i 

with Li and A the length and breadth of the paral- 
lelepiped, while A) is the edge-length of the cube. For 
the HS mixture we have 
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where 



6 = \ Pi D i> 6 = 7T J] ft A 2 ' 



(E.7) 
(E.8) 
: , (E.9) 

(E.10) 



and A (i = 1> 2) the particle diameters. 

To calculate the r parameter we only need expressions 
for en and ci 2 . The Fourier transforms of the geometric 
measures of overlapping bodies for the HPAR-HC fluid 
with a wave number k = (0, 0, k) (smectic symmetry) are 



/n(A) = 8A! 

Rn(k) = 4A 2 



2 sin(fcLi) 



k ' 

sin(fcLi) sin(fcA) 
u \ ; 1 U X ; 



k 

sin(ALi/2) 



(E.ll) 



(E.12) 



S n (k) = 2 A 3 



, sin(A:Li/2)\ 2 ^ fsm(kL 1 /2)' s 2 



D 



sin(fcLi) 



Vn(g) = 4 A 4 



4 / sin(AA/2) 



(E.13) 
(E.14) 
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The expressions for particles 1 and 2 are 
Mk) 

R 12 {k) = 2D 12 



2 sin(fcLi 2 /2) 
zlJ i2 r 



(E.15) 
sin(fcLi 2 /2) 



2Di 



sin(/cLi/2)sin(fcD 2 /2) 



Sia(g) = 2AA 



2D 



12" 



A- 2 



sin(fcLi/2)sin(fcD 2 /2) 



(E.16) 



2A 2 
AA 



sin(/cLi/2)sin(/cL»2/2) 



fc 2 

sin(fcLi 2 /2) 



V 12 (q) = 4(AA 



2 sin(fc J Li/2)sin(fcD 2 /2) 



(E.17) 
(E.18) 



where D\ 2 = A + D 2 and A2 = L\ + A- 

For HS mixtures these expressions have the form 



fijik) 



4ir 

fc 2 

4tt 
fc 2 



sin(fcAj/2) £> 



cos(*Aj/2) 



(E.19) 



sin(fcA/2)sin(fcD i /2) 



^cos(fcAi/2) 



(E.20) 



SVj (fc) — 



fc 2 



Aj sin(fcA/2) sin(fcA/ 2 ) 



A: 2 



AA sin(jfeAj/2) 



(E.21) 



4 fc 

(4tt) 2 rsin(fcA/2)sin(fcA/ 2 ) sin(fcAj/2) 



fc 4 



A; 2 



(E.22) 



+ cos(fcA/2)cos(fcA7 2 )] i 
for i,j = 1,2. 

Appendix F. Asymptotic behaviour of r 



The asymptotic behaviour of (dA/dp 2 )^ for = 
A/A > 1 can be calculated from ()C.8|) and (jE.lj) . 
particularizing the latter equation to HPAR-HC and HS 
mixtures. After some algebra, we arrive at 



<9A 

dp 2 



(0) 



9cn 

9P2 



(0) 



2 / HPAR 

^y(l + y){2(l + y) 2 T?(k* /2) 



[T (k*) + 2yT 1 (k* ) + 2/ 2 T 2 (fc*/2)] nf) (F.l) 



where Ki = A/ A, ?o(x) = cosx, and T\(x) = sinx/x, 
while y = rjo/(l — rjo), and fcp = fcoA- The expression 
corresponding to the HS fluid reads 

(0) 

<Wffs 



^ ' - -^y(l + y){2y 2 (l + 3y)T v (k* ) + T s (k*) 



2yT f (k* ) + y(l + 6y)T s (k*)}. 



(F.2) 



where we have defined the dimcnsionlcss quantities 
T v {kl) = 6v^ 2 V(k* ), T s (k* ) = 3a?6(k$), T f (k* ) = 
3u 1 1 /(fco), and Ts(/cq) = 6(wiSi) _1 S'(A:q), with si and ui 
the surface area and volume of particle 1, while S(q) is the 
Fourier transform of the Dirac delta function (5(A — |r|). 
The presence of a delta function indicates that, in the 
limit d — + oo, the one-component sticky-sphere limit of 
the fluid is obtained. We arrive at the same conclusion, in 
the same limit, for the HPAR/HC mixture, by compar- 
ing Eqns. (jF.ip and (|F.2[) . Each term in the right-hand 
side of Eqn. (jF.ip . from left to right, has the same mean- 
ing as in the HS mixture: they are related to the Fourier 
transform of the volume, adhesiveness, overlap function, 
and surface area of the overlap region between two parti- 
cles, respectively. For example, the inverse Fourier trans- 
form of T (fc*), with the smectic symmetry q = (0,0, q), 
results in a term proportional to <5(A — |z|)/2, which 
reflects the stickiness of parallelepipeds along the z di- 
rection. A complete effective density functional for the 
infinite asymmetric limit was worked out for hard-cube 
mixtures in Ref. [23j| . The direct correlation function re- 
sulting from this functional has a Dirac delta function 
located at contact of the sides of the cubes. 

However, there is an important difference between the 
HP-HC and HS mixtures, which is related to the square 
and cubic power dependence of the expressions (|F.1[) 
and (|F.2[) with respect to the asymmetric parameter 
d = A/ A- This result is related to the fact that the 
planar geometry of the sides of parallelepipeds enhances 
the depletion interaction between two big particles. We 
should compare these results with those obtained using 
the Onsager approach, where Cjj(fc, pi, p 2 ) — /y(fc) and 

then {dA/dp 2 ) {0) ~ O(d ). Thus, we can conclude that 
depletion interactions in this limit cannot be properly 
described by the Onsager approximation. 

Finally, to obtain the asymptotic behaviour of the r 
parameter for the HP-HC mixture in the limit d — > oo 
, we use Eqns. (|C.10p . (|F~T|) and (|EJ]) . particularized 
for HP, together with the spinodal instability condition 

0, to obtain explicitly 



1 - Poc\i 



KidW(ni), 



(F.3) 



W{k x ) = - {2(1 + y) 2 T?(k*/2) + [T (k*) + 2yT 1 (k%) 
+ y 2 T 2 (k*/2)} K r 2 } {y- 1 (3 + y- 1 ) + 2(3 + y^k*) 
_ [ 3 + y + 6(l + y) 3 ]r l 2 (fc*)}- 1 (F.4) 
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In Fig. JTTJ the function W(«i) is plotted. As can be FIG. 11: The function W(m). 

seen from the figure, (i) W varies very little with respect 
to Ki, and (ii) depletion is maximum in the Onsager limit 
(«i — > oo). Note that, on taking the latter limit, the 
condition d — D1/D2 ^> 1 must be fulfilled. 

It is well known that, when the PY approximation is 
used, the condition 1 — poCn(k, po) > is always fulfilled 
for all r\ and q in the physical parameter region. Because 
of the absence of a fluid-solid spinodal in this approxi- 
mation, we cannot calculate the r parameter to estimate 
the effect of depletion at the freezing transition. 
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